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We study how the spontaneous relaxation of a qubit affects a continuous quantum non-demolition 
measurement of the initial state of the qubit. Given some noisy measurement record "J, we seek 
an estimate of whether the qubit was initially in the ground or excited state. We investigate four 
different measurement protocols, three of which use a linear filter (with different weighting factors) 
and a fourth which uses a full non-linear filter that gives the theoretically optimal estimate of the 
initial state of the qubit. We find that relaxation of the qubit at rate 1/Ti strongly influences the 
fidelity of any measurement protocol. To avoid errors due to this decay, the measurement must be 
completed in a time that decrease linearly with the desired fidelity while maintaining an adequate 
signal to noise ratio. We find that for the non-linear filter the predicted fidelity, as expected, is 
always better than the linear filters and that the fidelity is a monotone increasing function of the 
measurement time. For example, to achieve a fidelity of 90%, the box car linear filter requires a 
signal to noise ratio of ~ 30 in a time T\ whereas the non-linear filter only requires a signal to noise 
ratio of ~ 18. 

PACS numbers: 03.65. Yz, 03.67.Lx, 03.65.Wj, 89.70,+c 



I. INTRODUCTION 

In this paper we consider the following problem. Given 
a qubit initially either in the ground or excited state 
with finite lifetime Ti, how can we best make use 
of a continuous-in-time noisy quantum non-demolition 
(QND) measurement to optimally estimate the initial 
state (ground or excited) of the qubit? A number of 
authors have previously considered the related prob- 
lem of optimal estimation of the present state of the 
qubit based on the past and current measurement record 
P, [H, H, HL S but, to the best of our knowledge, the 
problem we consider here has not been previously stud- 
ied. 

QND measurements play a central role in the theory 
and practical implementation of quantum measurements 
Q. In a QND measurement, the interaction term in the 
Hamiltonian coupling the system to the measuring ap- 
paratus commutes with the quantity being measured, so 
that this quantity is a constant of the motion. This does 
not imply that the quantum state of the system is totally 
unaffected, but it does imply that the measurement is re- 
peatable. For example, a Stern-Gerlach measurement of 
a z for a spin 1/2 particle initially prepared in an eigen- 
state of a x will randomly yield the results +1 and —1 
with equal probability. However, all subsequent mea- 
surements of <7 2 will yield exactly the same result as the 
initial measurement. 

The fact that QND measurements are repeatablc is of 
fundamental practical importance in overcoming detec- 
tor inefficiencies. A prototypical example is the electron- 
shelving technique [1, Q used to measure trapped ions. 
A related technique is used in present implementations 
of ion-trap based quantum computation. Here the (ex- 
tremely long-lived) hypcrfine state of an ion is read out 



via state-dependent optical fluorescence. With properly 
chosen circular polarization of the exciting laser, only 
one hyperfine state fluoresces and the transition is cy- 
cling; that is, after fluorescence the ion almost always 
returns to the same state it was in prior to absorbing 
the exciting photon. Hence the measurement is QND. 
Typical experimental parameters [Toj allow the cycling 
transition to produce N ~ 10 6 fluorescence photons. 
Given the photomultiplier quantum efficiency and typi- 
cally small solid angle coverage, only a very small number 
fid will be detected on average. The probability of getting 
zero detections (ignoring dark counts for simplicity) and 
hence misidentifying the hyperfine state is P(0) = e~ nd . 
Even for a very poor overall detection efficiency of only 
10~ 5 , we still have fid = 10 and nearly perfect fidelity 
F = 1 - P(0) ~ 0.999955. It is important to note that 
the total time available for measurement is not limited 
by the phase coherence time (T 2 ) of the qubit or by the 
measurement-induced dephasing [H, H, [TJ [l2| , but rather 
only by the rate at which the qubit makes real transi- 
tions between measurement (er z ) eigenstates. In a per- 
fect QND measurement there is no measurement-induced 
state mixing Q and the relaxation rate 1 jT\ is unaffected 
by the measurement process. 

The ability to read out a qubit with high fidelity is 
of central importance to the successful construction of a 
quantum computer [l^ . In order to successfully mea- 
sure a qubit, its quantum state must be mapped into a 
piece of classical information by measuring the relative 
occupation of its two states with the highest possible fi- 
delity. Possible qubit implementations include supercon- 
ducting circuits, silicon based electron and nuclear spins, 
and trapped ions, among others [3, [HI, [l|| Qj], HE E, 
I20L l2ll , I22I [23j . In order for qubits prepared in different 
states to be distinguishable, the measurement must be 
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completed before the excited qubits decay [J, [24|. Many 
atomic qubits have sufficiently long lifetimes so that re- 
laxation is not a major concern [HI, EH, HO], but most solid 
state qubits have lifetimes on the order of microseconds 
or less, and spontaneous relaxation plays a significant role 
in the measurement. The qubit relaxation affects differ- 
ent measurement schemes differently, but in all cases, it 
can limit the maximum fidelity. 

Although the behavior of a qubit during continuous 
measurement has been studied using Monte Carlo simu- 
lations @, S I3j IB- 6j, no attempt has been made to derive 
an analytical expression for the probability distribution 
of the initial state or to study how spontaneous emission 
impacts the measurement fidelity. This is at least par- 
tially because there exist very few high fidelity contin- 
uous QND measurements, and even fewer that operate 
in a regime where qubit relaxation is a limiting factor. 
With the application of low temperature amplifiers, high 
fidelity (though not necessarily weak continuous) mea- 
surements of superconductin g qu bits are now becoming 
feasible [H, M, HE M, Hi HE M, ED, E2] . These mea- 
surements have found asymmetries between the probabil- 
ity distributions for the integrated signal corresponding 
to the ground and excited states, but could not accu- 
rately predict them. 

Here we find that when we use a measurement pro- 
tocol that only records the integrated signal, the qubit 
relaxation induces asymmetry in the probability distri- 
butions, and that with a sufficiently precise detector, 
the distributions become distinctly non-Gaussian. Un- 
like measurement of a perfect (i.e., non-decaying) qubit, 
where fidelity is always improved by a longer measure- 
ment, we show that there is some optimal measurement 
fidelity that depends on the signal to noise ratio (SNR) 
of the detector and the filter used. The first filter we 
consider is the linear box car filter and optimize over the 
integrated time tf. Choosing a longer or shorter inte- 
gration time will lower the fidelity of the measurement. 
Next we show that by choosing a filter that gives expo- 
nentially less importance to results at later times slightly 
increases the fidelity. We then numerically find the opti- 
mal linear filter and compare these linear filters to a non- 
linear filter that yields the theoretically optimal estimate 
of the initial state of the qubit given some measurement 
record 'J. We find that we can reach the same fidelity as 
the linear filters at a substantially lower SNR. Further- 
more, due to the nature of the updating protocol, the 
fidelity is a non-decreasing function of the measurement 
time. In summary, in this paper we determine the opti- 
mal measurement fidelity given four measurement proto- 
cols for continuous measurement experiments currently 
being performed, and also provides a guideline for the 
necessary detector signal to noise ratio in order to reach 
a particular desired fidelity in future experiments. 

There are two major ways of measuring qubits. The 
first method is a latching measurement, for example by 
having the qubit state modify the switching current (or 
state) of an adjacent Josephson junction [28|, [29[ or the 



bifurcati on p oint of the non-linear Josephson plasma os- 
cillation [3d [31 , 32]. In such latching measurements, the 
qubit is measured very quickly with very hi gh s ignal to 
noise, and after only a short waiting time (l5l Il7l.l28l.l32l]. 
Some versions of such strong measurements can in prin- 
ciple be QND [3^. The second method is to perform 
a sequence of repeated or weak continuous quantum 
non-demolition (QND) measurements, which each leave 
the populations of the qubit unchanged. Several recent 
experiments with solid-state qubits [23|, [H, Hj], have 
used continuous QND measurement schemes in which the 
qubit lifetime imposes the main limitation on the mea- 
surement fidelity, for which the analysis of this paper 
should apply. 

It is not just continuous measurements that are af- 
fected by qubit relaxation. For example, consider an 
idealized latching measurement scheme where a qubit is 
prepared in an eigenstate, but there is some finite arming 
time t arm before a perfect measurement is made. In this 
case, a qubit prepared in the ground state will always be 
measured correctly, but a qubit prepared in the excited 
state may have decayed during i arm and be misidenti- 
fied. Thus if iarm is not infinitcsimally small, the qubit 
lifetime places a limit on fidelity even with a perfect de- 
tector. 



II. QUBIT WITH INFINITE LIFETIME 

It is worthwhile to first consider the case when the 
qubit cannot relax from the state it is initially prepared 
in, so it is "fixed" for all time. This allows us to formalize 
our intuitive understanding of a general continuous mea- 
surement and to give us a result against which we can 
compare the finite lifetime case. We consider the mea- 
surement of a qubit with two states |+) (excited) and 
|— ) (ground) and assume that the measurement result is 
given by the actual value of the qubit state plus Gaus- 
sian noise. This assumption is justified for example in 
the current circuit quantum electrodynamic experiments 
[HI H2I I22I , |23| in which a cavity is dispcrsivcly coupled 
to the qubit (no energy is exchanged between the cavity 
and the qubit). A homodyne measurement on the cavity 
output will reveal the cavity state (which is proportional 
to a z ) plus Gaussian noise [H|,[2l|]. This Gaussian noise 
will be at least the photon shot noise but in present ex- 
periments it is dominated by the following amplifier. In 
other words, the measurement is faithful and given that 
the system is in state i = ±1 our detector for a time 
interval dr outputs ?P(t) with statistics 



f (VI*) = \ „ exp[-(V>- ifdr SNR/2]. (2.1) 
V 2n 

For convenience we have also introduced a dimensionless 
time r = t/T± where T\ is an arbitrary but finite number, 
that will become the relaxation lifetime when we treat the 
finite lifetime case. Here SNR is the ratio of integrated 
signal power to noise power. It is linear in the integration 
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time and we will adopt the convention of specifying the 
SNR as that achieved after integrating for time Ti. 

From this distribution we can write ip(r) in terms of 
the Wiener increment dW(r) [3{| as 

ip{T)dT = i±(r)dr + V SNR _1 dW(r). (2.2) 

Here we have introduced the subscript ± to indicate a 
possible realization of the dynamics of the qubit given 
the initial condition ±1. For this case the qubit can be 
initialized in either state, but because it has an infinite 
lifetime, it is fixed in whatever state it starts in for the 
duration of the measurement. That is, i±(r) = ±1. 

We define our measurement signal s as the output of 
the detector integrated over time Tf 

s= f drip(T). (2.3) 
Jo 

Formally we are restricting ourselves here to a simple box 
car linear filter which uniformly weights the measurement 
record in the interval < t < T[. Using Eq. (|2.2[) it 
is simple to carry out the above integral and rewrite the 
measurement signal as s±(t{) = ±Tf + Xq[0, a 2 ] where 
Xq[Q, a 2 ] is a Gaussian random variable of mean and 
standard deviation a = y^/SNR. The signals then fol- 
low the familiar Gaussian distributions 

u \2tt 

Because these distributions are symmetric about s = 0, 
the most obvious analysis is to set a signal threshold 
Vth = and call every measurement with s > v t h a (+1) 
state, and every measurement with s < vth a ( — 1) state. 
Calculation of fidelity in this case is accomplished using 
the definition 

dsP+(s) - / dsP-(s). (2.5) 

For the case of infinite qubit lifetime we have the simple 
result 

F-«f(V^). (2.6) 

A fidelity of zero corresponds to a completely random 
measurement that extracts no information, a fidelity of 
one corresponds to a perfect faithful measurement, and 
in between the measurement conveys varying degrees of 
certainty. As Tf becomes large, Eq. (|2.5[) predicts that the 
fidelity rapidly approaches unity. Higher SNR serves to 
speed up the convergence, but as long as SNR is non zero, 
any desired fidelity is attainable simply by measuring the 
qubit for long enough. In TablelH the required SNRfi xe( j is 
listed in order to achieve a given fidelity within T\ . Note 
the same results for the fidelity would be obtained if we 
used the optimal non-linear filter of Sec. IVII That is, for 
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FIG. 1: (Color online) Example detector output as a func- 
tion of time for a continuous QND measurement with SNR = 
570, (dr — 1 x 1CP 4 ) which allows for a measurement fidelity 
of 99%. The dashed red line corresponds to a large ensem- 
ble average of qubit outputs, and decays exponentially with 
characteristic time Ti . The solid blue line corresponds to the 
measurement record for a particular single shot. At this value 
of SNR, the jump in output is clearly visible and occurs at 

a qubit with infinite lifetime the simple box car filter is 
equivalent to the optimal non-linear filter, it is only when 
we include relaxation this is not the case. This will be 
discussed in detail in the next four sections. 



III. BOX CAR LINEAR FILTER 
A. Probability distributions for QND measurement 

Here we consider the same measurement protocol as 
in Sec. [TTJ however unlike a qubit with infinite lifetime, 
where both states behave similarly, a qubit with a finite 
lifetime has a fundamental asymmetry in how the excited 
and ground states behave. We assume (in our model) 
that the excited state decays at rate 1 jT\ but the transi- 
tion rate upward out of the ground state is zero. A qubit 
prepared in the ground state will never experience any 
excitations, so it can be treated much the same as the 
fixed qubit discussed above and hence P-{s) = P^ xed (s). 
By contrast, an initially excited qubit will produce an 
ensemble averaged output which will decay exponentially 
with characteristic dimensionless time ti = 1. Although 
most qubits dephase in some shorter dimensionless time 
T2, our only concern here is the relative population of the 
two qubit states, so we are not limited by decoherence of 
the qubit. For a single qubit, this translates into a sin- 
gle, abrupt relaxation of the qubit at some dimensionless 
time Ta that is exponentially distributed with mean di- 
mensionless time 1, P{Td) = exp(— Td). That is, if the 
qubit was initially in the excited state then i±(r) would 
obey 

*±(T)=0(T d -r)-0(T-T d ). (3.1) 
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Thus given a possible realization i±(r), we can generate a 
typical record an experimentalist would measure by using 
Eq. ()2.2j) . A typical trajectory for tp{r) is shown in Fig. [I] 
for a SNR of 570. Here we see, at this value of SNR, the 
jump in output is clearly visible and occurs at ^ 1 

In the case where the qubit happens to decay early, 
Td <C 1, the signal s from a qubit initially prepared in 
the excited state would be almost indistinguishable from 
that of a qubit initially prepared in the ground state, and 
even if the measurement apparatus were nearly perfect, 
almost no information could be extracted. The proba- 
bility distribution P+(s) for the measurement signal of 
an initially excited qubit can be determined analytically 
with a simple derivation. This will allow a more quan- 
titative discussion of fidelity, and will eventually allow 
for optimization of integration time, Tf. The critical dif- 
ference between the infinite lifetime system described in 
Sec. Inland an actual QND measurement is that s is now 
a function not only of the dimensionless integration time 
Tf, but also the exponentially distributed random relax- 
ation time Td. 

If the qubit is initially in the excited state, then from 
Eqs. (|2.2[) . (|2.3p . and (|3. If) a possible realization of s± 
will be 

S± = T[6(T d - Tf) + (2t c1 - T{)9(n - Td) 

+X G [0,a 2 ]. (3.2) 

That is, the probability distribution for s given a realiza- 
tion with a decay at time Ta is 

P+(s\r d ) = -^exp[-(.s-T f )72a 2 MT d -Tf) 
V27rer 

+ -^cxp{-[ S -(2T d -T f )] 2 /2a 2 } 
\>2ir<r 

x6(T{-r d ), (3.3) 

and from this one can easily obtain the probability dis- 
tributions for s by averaging over all possible realizations 
(decay times). Doing this gives 

P + (s) = c -(s-T f ) 2 /(2a 2 )-r f + lg-(«+T,)/2 

><^ /8 H^^)-erf(^^)} (3.4) 

Although it does not figure directly into this analysis, 
it is not difficult to expand this treatment to consider 
a measurement with finite "demolition" that stimulates 
both excitation and relaxation of the qubit. In this case, 
both P_|_ and P_ will be non-Gaussian, because the qubit 
may excite and relax several times during the measure- 
ment interval. To calculate the distributions, all we need 
to do is extend the possible realization to include multi- 
ple relaxations and excitations keeping in mind that the 
relaxation and excitation times are not independent vari- 
ables: Td n must occur before T dn+1 . For stronger or less 
ideal measurements, it may become necessary to include 
these extra terms, but for now the demolition is taken 
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FIG. 2: (Color online) Evolution in time of probability dis- 
tributions for QND measurement of qubits initialized in the 
ground (blue) and excited (red) states for the box car linear 
filter and a SNR = 30. For short Tf, the distributions are 
almost symmetrical because an excited qubit has probably 
not decayed yet. As Tf gets longer, the qubit is much more 
likely to have decayed, and the mean of P+(s) begins to drift 
back towards that of P-(s). At some optimal point in be- 
tween, Tf = T opt , and if the threshold i/fh is chosen so that 
P+(iAh) = P-(vth), the misidentified tails are minimized, and 
fidelity is maximized. 

to be small compared to the spontaneous relaxation and 
can be ignored. 

The Gaussian first term in Eq. (|3.4|) is dominant for 
Tf < 1, and P+(s) is nearly symmetrical to P_(s). For 
these fast measurements, there is little chance that the 
qubit decays during the measurement, and the distribu- 
tions will be very similar to those of fixed state bits. At 
some point, the probability that the qubit has relaxed 
will be large enough that it significantly affects the dis- 
tributions. At this point, predictions based on the as- 
sumption of no relaxation are no longer valid, and the 
non-Gaussian term in Eq. (|3.4[) becomes very important. 
This effect can be seen in the non-Gaussian tail of P+(s) 
in the second two time cuts in Fig. [2] 

A strongly non-Gaussian P + distribution is not re- 
quired for the fidelity to exhibit a maximum at some 
finite measurement time. Even though for experiments 
with SNR w 1, the non-Gaussian tail of P + is not very 
prominent, the asymmetry between the distributions in 
both height and width is quite clear. For long enough 
dimensionless integration times, a qubit initially in the 
excited state will have probably relaxed relatively early 
in the measurement, so its mean should be very similar 
to that of a qubit initially in the ground state, the distri- 
butions will be almost identical, and all resolving power 
will be lost. 



B. Optimal box car filter 

The behavior of fidelity as a function of integration 
time for a QND qubit measurement is very different from 
a measurement of a fixed state qubit that never relaxes. 
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Recall that in the fixed state case, the fidelity eventu- 
ally converges to one, independent of SNR. We define fi- 
delity as we did for the fixed state case in Eq. (|2.5p . The 
difference here is that because the distributions are not 
necessarily symmetrical, the signal threshold vth is not 
always zero. Maximizing F with respect to vth yields the 
following implicit equation for Vth 



P+(w th )=P_(i^ h ). 



(3.5) 



As an aside we note that from this we see that we can 
write the fidelity (optimized with respect to vth) in the 
alternative form 



F 



ds\P + (s)-P-(s)\ 



(3.6) 



Despite the complications of finite lifetime, the integra- 
tions in Eq. (|2.5p can still be carried out analytically to 
yield 



F 



e - 2 /8 e -(H> 



l 1 +r)/2| 



erf 



-erf 



a 2 - 2(;/ lh + r) 
2V2a 



a 2 - 2(i/ th - t] 
2V2a 

} (3-7) 



The fidelity is maximized by numerically solving for u t h 
such that P_(^th) = P+{vt\\)- Using the correct value of 
^th(Tf); it is straightforward to compute F(j{) and then 
vary Tf to obtain the optimal value of the integration 
time. 

Roughly speaking, fidelity is a measure of how separate 
the probability distributions are, ranging between for 
complete overlap and 1 for no overlap. Fidelity is limited 
by detector noise for short times and spontaneous relax- 
ation for long times. For Tf < 1, the probability of a 
relaxation during the measurement is very low, and the 
fidelity behaves similarly to the fixed state case: it in- 
creases with increasing Tf. For longer 7? »1, the qubit is 
more and more likely to have relaxed during the measure- 
ment and the mean of P+(s) will stop increasing linearly 
and in the long time limit will actually decrease in time 
as can be inferred from the decrease in optimal threshold 
value plotted in Fig. [3]B) (red solid line). This implies 
the existence of some intermediate time T opt (SNR) that 
maximizes the fidelity, this is clearly seen in the red solid 
line of Fig. [3] A) where the fidelity for a SNR= 10 has a 
maximum of 0.79 at T opt = 0.34. 

To find this time for a given value of SNR, we com- 
pute F(t{), then numerically solve dF(T{)/dri = for 
Tf = Topt- The solution, F opt = F(T opt ) is the best 
possible fidelity for the simple linear box car filter; no 
improvement can be made from this value without im- 
proving the measurement apparatus or the qubit lifetime. 
This optimal fidelity is only achievable by correctly set- 
ting Tf = T opt and v t h. Any variation of these parameters 
will reduce the fidelity. It should be pointed out that 
while the condition P-{v t \ x ) ~ P + (is t h) does maximize 
the fidelity, it implies that the measurement protocol is 
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FIG. 3: (Color online) A) The fidelity F as a function of 
time Tf for all measurement protocols. The box car filtered 
integrated signal is the red (solid) line and the exponential fil- 
tered integrated signal is the blue (dashed) line. Here we see 
that both these schemes have an optimal measurement time 
and that the latter case is less sensitive to the measurement 
time. The non-linear filter is shown by the green (dashed- 
dotted) line. Here we see that it is clearly better than the 
other cases and that there is no optimal measurement time. 
B) The threshold v t ^ as a function of time Tf for both the box 
car filtered integrated signal (red solid line) and the expo- 
nential filtered integrated signal (blue dashed line). In both 
subplots the SNR after time Tf = 1 is 10 and time is measured 
in units of Ti . 



biased towards the ground state. That is, if the qubit is 
prepared in the ground state we are more likely to assign 
it correctly than if it was prepared in the excited state. 
An unbiased measurement protocol is obtained by set- 
ting j/ t h such that ds P-{s) = f v ds P+(s). Doing 
this results in a slightly lower optimal fidelity. 

Since F op t depends only on SNR, it is possible to first 
derive the required signal to noise ratio after one lifetime, 
and then choose the correct measurement time and signal 
threshold required in order to attain any arbitrary fidelity 
as shown in Fig. [4] (red solid line in part A) along with 
the optimal measurement time (red solid line in part B). 
For example, one standard initial goal is a fidelity of 90%, 
sufficient to violate Bell's inequalities (36J. As shown in 
table HI this fidelity requires a minimum SNR of 30 after 
time T\. 

The following argument shows how rapidly the re- 
quired SNR diverges for very high values of fidelity. For 
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TABLE I: The required minimum signal to noise ratio SNR after T\ and the optimal optimal measurement time, T opt for a 
continuous measurement of a qubit with infinite lifetime (SNRg xc d), a simple box car linear filter (SNRbc), an exponentially 
decaying linear filter (SNR cxp ), the optimal linear filter (SNRol), and a non-linear Bayesian filter (SNRnl). The last column 
is the maximum allowed waiting time r arm for an idealized latching (instantaneous) measurement with infinite signal to noise 
performed after time r arm in order to achieve the desired fidelity F. 
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where x = In ( 1+ ^ R j . 

Approximating the error function in Eq. ()3.8[) by unity 
and neglecting ln(xo) relative to xq leads to the following 
simple asymptotic form for large SNR 
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(3.10) 



This can be rewritten as in terms of the optimal mea- 
surement time as Fq ~ 1 — T opt /2. Here we see that 
to achieve the optimal, the measurement must be com- 
pleted in a time that decreases linearly with the desired 
fidelity. The slight deviations form this zeroth order re- 
sult is shown in the inset of Fig. [4] B. This result is 
consistent with the intuitive picture that to achieve a fi- 
delity F = 1 — rj with ?/ <C 1, we must have sufficient SNR 
to be fooled only by those extremely early decays which 
occur in time r* ~ 2rj. This set of decays occurs with 
probability ~ 2ry and is the main source of the infidelity. 



FIG. 4: (Color online) A) Optimal fidelity F opt and B) mea- 
surement time r op t as functions of SNR. The inset in plot 
B shows the optimal fidelity as a function of the measure- 
ment time. The red (solid) line is for the box car filter, blue 
(dashed) is for the exponential filter, black (dotted) is for the 
optimal linear filter, and the green (dashed-dotted) is for the 
non-linear filter. Here we see that by using the optimal non- 
linear filter we get an improved fidelity in comparison to the 
simple linear filters. Particular values are shown in Table [T] 



large SNR, it is a good approximation to set the thresh- 
old f t h = and then the fidelity is approximately 

F = cxp(-T f /2)crf (ySNRrf/2) (3.8) 

Optimizing this with respect to Tf and then using the 
asymptotic form for the error function of large argument 
yields the following expression for the optimal integration 



IV. EXPONENTIAL DECAYING LINEAR 
FILTER 

In this section we consider the case when we include 
an exponential weighting factor in our integrated signal 
s. This is chosen as this would be the optimal linear fil- 
ter if our signal was simply a decaying exponential with 
a random initial amplitude, this can be proven by min- 
imizing our estimate of the amplitude in a least square 
sense. However, even though our signal on average is of 
this form, in one particular run it is not. Thus this will 
not be the optimal linear filter for estimating the initial 
state of the qubit. For this reason we introduce an in- 
tegration time Tf and optimize the measurement fidelity 
over this. That is, the signal s for this filter is given by, 

s = [ dr 7/>(T)e~ T . (4.1) 
Jo 
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As in the last section we first determine a possible s given 
that the qubit started in the excited state. Doing this 
gives 



VSNR 



Jd - t) 

■Tf 



dW{r) 



(4.2) 



o 



and by treating the noise integral as simply a linear com- 
bination of infinitesimal Gaussian variables gives 



1 - 2e 



~Td , „-T{ 



i T l - T d) 



+a 0(r d -r f )+X G [O, ( 7 2 ], 



(4.3) 



where a = - e- 2r f)/2 SNR and a = 1 - e~ Tf . Here 
we see that unlike before as the measurement time be- 
comes large the variance saturates at (2 SNR) -1 rather 
then continuing to increase linearly with time. That is, 
we have designed our filter such that when all the infor- 
mation about the qubit state has been lost into the T\ 



J 



environment the noise in the integrated signal will remain 
constant. 

Following the same procedure as before the excited 
state distribution is 



P+(s) 



erf 



f a + s 
V V2a 



erf 



y/2o 



2ira 



exp 



(g - a Y 
2o 2 



n 



(4.4) 



Repeating the above for the ground state initial condition 
gives st = — o,+Xg[0, a 2 ] and a ground state distribution 
of the form 



P-(s) 



2n<r 



exp 



(s + a)" 
2a 2 



(4.5) 



With the above two distributions the fidelity as before, 
can be solved analytically and in terms of the v t h is 



F 



exp 



(a - v th f 



2a 2 



exp 



(a + v t h) 
2<r 2 



(e- 



r 



1 - Vth) 



erf 



a - vth 
V2a 



erf 



a + Vth 
V2cr 



(4.6) 



Using the same procedure as before we can numeri- 
cally determine the vth which maximizes the fidelity For 
a SNR of 10 (in time T\) the fidelity as a function of mea- 
surement time is shown in Fig. [3]A) as a blue dashed line. 
Here we see as before there is an optimal measurement 
time. To measure any longer than this time results in 
a lower fidelity. The optimal fidelity and measurement 
time are shown in Fig. 0] (blue dashed line) as a function 
of the SNR. Here we see that by using this filter, the fi- 
delity is slightly better than the simplest case (see TableU 
for some values), but more importantly the curvature of 
the fidelity at r opt is less. This means that this filter is 
less sensitive to errors in the measurement time. That is, 
this protocol would be more practical to implement than 
the simple box car integrated signal of Sec. IIIII 



is defined as the difference between the probability of us 
making a correct assignment and an incorrect assignment 
[Eq. (|2.5|) ]. The assignment criteria we use is again if s is 
above z/ t h then we say the qubit was initially up and if it 
is below vth then it was down. Vth like k(t) is determined 
by the maximization procedure which we will describe 
now. 

For an unknown kernel the signal conditioned on the 
qubit being initially up will be given by 



,s+ = 2a(r d ) -a(oo) +X G [0,a 2 }, 



where 



o(I) 



fc(r')dr' 



V. OPTIMAL LINEAR FILTER 



k 2 (T')dT' /SNR. 



(5.2) 

(5.3) 
(5.4) 



In this section we calculate the optimal linear filter for 
estimating the initial state of the qubit. We define the 
linear signal s by the relation 



k(r)ip(r)dT, 



(5.1) 



where the kernel k(t) is unknown and is determined by 
maximizing the measurement fidelity. This, as before, 



If the qubit was initially in the ground state then the 
signal would be 



s_ = -a(oo) +X G [0, a 1 



(5.5) 



^From the above equations the excited state and ground 
state probability distribution for the signal s are 



P+{s) = f 
Jo 



cxp[-r d - (s - 2a(r d ) + a(oo)) 2 /2cr 2 ] 

\l2"K0 2 



dTA 
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P-(s) 



exp[-(s + a(oo)) 2 /2cr 2 



V27TCT 2 

and by using Eq. I|2.5p the fidelity is 



(5.6) 
(5.7) 



F = - 



f e "" erf ( 



2a(r d ) - a(oo) - ^ th 



\72o- 



dr d , 



1 / -a(oo) - z^th ' 

2 CT V V2a 



(5.8) 



Maximizing this gives the following set of coupled differ- 
ential equations 



d T a(r) 
d T k(r) 



fc(r), 
-exp{- 



2a(r)[a(r) — a(oo) 



(5.9) 
t'thl/o' 2 }, 
(5.10) 



with the initial conditions A:(0) = 1, a(0) = 0, and the 
boundary condition k(oo) = 0. ft is this latter condition 
which determines i^th- This system of equations can be 
solved numerically using a shooting method [33] ■ The re- 
sults arc shown in Fig. [5] for a SNR of 1.0 and 172. Here 
we see that for the small SNR, the kernel k{r) can be ap- 
proximated well by exp(— (3t) where (3 is a fit parameter 
that is approximately equal to 1 + SNR/2. In the large 
SNR limit the k{r) cannot be fit by an exponential. For 
illustrative purposes the optimal linear filter is compared 
with the optimal box car linear filter in Fig [5] B) for a 
SNR of 172. Here we see that the time when the box car 
linear filter turns off, is comparable with the time scale 
of the optimal linear filter. 

Numerically solving for k(r) for a given a we can use 
Eqs. (|5.4[) and (|5.8p to plot the fidelity of the optimal 
linear filter as a function of the SNR. This is shown in 
Fig. |4] as a dotted black line and in Table [J as column 4. 
Here we see that the optimal linear filter out performs the 
other linear filters and is almost as good as the optimal 
filter (which is non-linear) and is described in the next 
section. 



VI. OPTIMAL NON-LINEAR FILTER 

In the previous sections we considered only the case 
where one looks at or only has access to the integrated 
signal, s(r). Here we assume that we have access to the 
full record W = {ip(r)} and ask how much better can we 
do with a non-linear filter. In particular, given this record 
what is our best guess at the initial state of the qubit. 
Mathematically, our best guess can be represented by the 
probability distribution P(io\^f). This is the probability 
that the initial condition io is ±1 given the record, \&. As 
with all probability distributions, this will range from 
to 1, and the closer it is to one the more we are certain 
that the qubit was in the initial state io = ±1. This by 
definition is the optimal protocol as it is the best estimate 



^-4 



O 




0.04 

Time, x 



FIG. 5: (Color online) The optimal linear filter for a SNR of 
1.0 blue (dashed) and SNR of 172 red (solid) as a function of 
time. The green (dashed-dotted) line in part B) is the optimal 
box car linear filter for a SNR of 172 (r opt = 0.046). 



of io given the complete set of information available. To 
find this distribution we use Bayes theorem, 



P(i |*) = 



pm )p(io) 

£ l0 P(*|io)P(<o)' 



(6.1) 



where P(io) is the initial probability distribution and 
P(^\io) is the probability that we would measure record 
^ given that the experiment was initially prepared in 
state io ■ We will assume that our experiment can prepare 
unbiased initial states and as such we take P(io ± 1) = 
1/2. The conditional distribution P(^\io) is the proba- 
bility that we would measure record, \E', given that the 
experiment was initially prepared in state io- Thus to 
calculate P(io\^) all we need to do is calculate P(^\io)- 
This is the non-trivial step, but can be determined by 



p(*|io) = ^p(*h' P )p(ip| 



(6.2) 



where p labels a possible realization of the qubit trajec- 
tory. That is, P(^|zo) can be determined by taking the 
ensemble average of P(^\i p ) for all possible realizations, 
i p ( i p is given by Eq. 
excited state or i p = - 
^From our simple model for the noise, Eq. (|2.1[) . P(^/\i p ) 
is simply a multiplication of many Gaussians each cen- 
tered on the instantaneous value of i v . 



if the qubit is initially in the 
1 if initially in the ground state). 
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For an excited state initial condition, the probability of getting record "J will be given by 
I 









/ dr d 




~ Td exp 


Jo 












e~ Tf exp 




- d 






Jo 



exp 



dr(ip{T) + l) 2 SNR/2 



dr{ifj{T) - l^SNR/2 



(6.3) 



Here the first term represents all possible trajectories that have had a decay in time Tf (the integration time), the 
second term represents all the trajectories that did not decay in this time, and A is the normalization constant. If 
the qubit were initially in the ground state then the probability of getting record \& will be given by 



P(*|i = -1) = Aexp 



- f ' c2t(V>(t) + l) 2 SNR/2 
Jo 



(6.4) 



Now that we have expressions for P(^\io) all that we need to do to get our best estimate of the initial state is to 
use Eq. (|6.1[) . Doing this gives 



P(i = l|¥) = 



1 



P{io 



Norm 
1 

Norm 



f ' dr d e- Td exp{ [s(r d , 0) - s(r, r d )]SNR} + e~ Tf exp[s(r f , 0)SNR] 
Jo 

exp[-s(r f ,0)SNR] 



where the two time integrated signal s(r, t') is 

s(t,t') = J T dr" ip(r") 

and the norm is simply 

Norm = / dr d e~ Td exp[(s(T d , 0) - s(r, r d ))SNR] + e" Tf exp[s(r f , 0)SNR] + exp[-s(r f , 0)SNR], 
Jo 



(6.5) 
(6.6) 

(6.7) 
(6.8) 



r 



Here we see that to solve this equation we need to 
evaluate a double integral over a stochastic process. This 
is impractical to solve numerically, however, as shown 
in the Appendix we can easily recast these integrals in 
terms of two sets of two coupled stochastic differential 
equations which require similar computational resources 
to that used with the linear filters. 

To show a typical trajectory for this estimated initial 
condition, we randomly generated records for both an 
excited and ground state initial condition. Rather then 
plotting both P(i = l|\l>) and P(i = we define 

z = P(i = 1|*) - P(i = -1|*) (this is the estimator 
that replaces s used in the linear filters), this will range 
from —1 to 1 and the closer it is to one of these limits, the 
more certain we are that the initial condition which gen- 
erated \& is this value. The results of this simulation are 
shown in Fig. [6] A). The red (solid) line corresponds to 
the case when the initial state was the ground state and 
the blue (dashed) line is for the excited state. Here we 
see that for this typical trajectory our estimate is fairly 



good at predicting the real initial condition and we are 
almost certain for the excited case. However this is only 
a typical trajectory and to get a better understanding of 
the predictability of this method we use the same fidelity 
measure as before, that is we subtract our wrong guesses 
from our correct guesses. To be more specific, we define 
the following assignment procedure: if z > then we as- 
sign the qubit as up and if z < we assign it as down. 
If z is equal to zero we ignore the result (or flip an unbi- 
ased coin to make our decision). Given this assignment 
criteria we can define the fidelity as 



lim 

M— >oo 



E - E + E ■ 

2+l>0 2+l<0 Z_i<0 



5_i >0 



(6.9) 



where M is the number of randomly generated "J for both 
+1 and —1 initial conditions. This is shown in Figs. OA) 
and [4] A) a s the green (dashed-dotted) line for M = 10 4 
(specific optimal fidelities are listed in Table H}. Here we 
see that the fidelity is always better than the other cases 
and that there is no optimal measurement time. That is, 
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0.5 1.0 1.5 

Time, x f [Ti] 




-1.0-0.8-0.6-0.4-0.2 0.0 0.2 0.4 0.6 0.8 1.0 

Estimate state, z(x f ^ oo) 

FIG. 6: A) A typical trajectory for the estimated initial 
condition, z(jf), given a record 9 which is randomly generated 
for an excited state initial condition (blue dashed line) and a 
ground state initial condition (red solid line) . B) A histogram 
of z(ji — > oo) for 10 4 trajectories for both an excited state 
initial condition (dark blue bars column 1) and a ground state 
initial condition (light red bars column 2). The SNR after 
T f = 1 is 10 and time is measured in units of T\ . 



other measurement schemes. For comparison, consider 
a latching measurement where the qubit state triggers a 
classical switching event in the measurement apparatus 
[28l |29| . [30L I31I . |32 |. Such a measurement has the ad- 
vantage that the detector state stays latched for a very 
long time, so that noise from subsequent amplification 
stages is completely negligible and the SNR is effectively 
infinite. We can roughly model this process as an instan- 
taneous measurement with no errors but a finite arming 
time t aim needed to set up the pre-latched state of the 
detector. We assume that the arming stage occurs after 
the qubit is prepared but before it is measured. The mea- 
surement of a qubit in the excited state will be wrong if 
the qubit relaxes before the measurement is made. The 
probability that the qubit relaxes during t arm rises ex- 
ponentially towards 1, so the fidelity falls exponentially, 
F(T arm ) = e~ Tarm . As shown in Table HI r arm must be 
slightly smaller than r opt . Although these times are dif- 
ferent in that r arm is a maximum value whereas r opt is an 
optimal value for a given SNR, they provide a valuable 
comparison between the measurement schemes. 

Of course, no latching measurement is truly instanta- 
neous with perfect accuracy in translating the qubit state 
into latching events. Moreover, a latching measurement 
usually induces some mixing of states, so the actual fi- 
delity may be lower than for an equivalent continuous 
QND measurement. Both measurement schemes have 
their own strengths and weaknesses, but in cither case, 
qubit decay can be a significant limiting factor on the 
fidelity. 



VIII. CONCLUSION 



unlike the box car filter and the exponential linear filter, 
the fidelity is a non-decreasing function of the integration 
time. 

It should be noted that as in the previous protocols, 
this is a biased measurement which favors ground state 
preparations. This can be seen by looking closely at Fig. 
[6]B). This figure shows a histogram of 10 4 estimates when 
the system is prepared in the excited state (dark blue 
bars, column 1) and the ground state (light red bars, 
column 2) for a SNR = 10. If we use our assignment pro- 
cedure and subtract the wrong guesses from the correct 
guesses for each prepared initial state we get a fidelity of 
0.92 for the ground state and 0.76 for the excited state 
(an average fidelity of 0.84). This is because when the 
qubit is prepared in the excited state, rare, early decays 
will fool the detector. 



VII. COMPARISON WITH LATCHING 
MEASUREMENT 

The challenges from qubit relaxation to attaining 
high fidelity in a continuous measurement also exist in 



We have examined the effect of qubit relaxation on 
the estimation of the qubit initial state (excited or 
ground) using a continuous-in-time noisy quantum non- 
demolition measurement for four different measurement 
protocols. In these protocols the measurement results 
are integrated with a box car linear filter, an exponen- 
tially decaying filter, an optimal linear filter and a non- 
linear Bayesian filter which by definition is the optimal 
theoretical filter. We found that in all these protocols 
there exists a theoretical limit on the measurement fi- 
delity. The determining factor of this limit is the signal 
to noise ratio of the measurement. Our results are sum- 
marized in Table U where we see that the non- linear filter 
reaches the same fidelity as the linear filters even for sub- 
stantially lower required signal to noise ratio. Lastly we 
compare the continuous quantum non-demolition results 
with latching measurements and found that there is a 
quantitatively different but qualitatively similar limit on 
the fidelity of latching measurements also due to relax- 
ation. The signal to noise ratio required to do successful 
qubit single shot quantum non-demolition measurements 
should be attainable in the near future, and there is no 
fundamental reason why significantly higher fidelity mea- 
surements cannot be performed. 
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APPENDIX A: NUMERICAL PROCEDURE 
USED BY THE OPTIMAL NON-LINEAR FILTER 



simulating Eqs. (|6.5|) and (|6.6|) which contains a double 
integral over a stochastic process This is not practi- 

cal numerically and a much better method can be imple- 
mented by deriving a set of coupled stochastic differential 
equations. To do this we start by rewriting Eqs. (|6.5[) and 
MM as 



P(io = 1|*) 
P(*o = -l|tf) 



P-i,i/> + P+i,ip 

P-l,4> 
P-i,i/> + P+i,ip 



(Al) 
(A2) 



In this appendix we present the method used to sim- 
ulate the optimal non-linear filter. This filter requires 



where 



P +1 j = e~ Tt exp[ S (r f , 0)SNR] + / dr d e~ rd exp{[s(r d , 0) - s(r f , r d )]SNR}, (A3) 

Jo 

P-W = exp[-s(r f) 0)SNR], (A4) 
and the signal s(t, t 1 ) is defined in Eq. (|6.7p . We first differentiate P+i,i/> with respect to Tf and by defining 

A+i.v, = e~ T{ exp[s(r f ,0)SNR] - ( ' dr d e~ Td cxp{[s(r d , 0) - s(r f , T d )]SNR}, (A5) 

Jo 



we get the following set of coupled stochastic differential 
equations 



— P+i^ = SNR V(t)A+i, v „ 

<2Tf 



(A6) 



d,T{ 



A+1.V, = SNRV(t)P+i^-(A+i^+P + i^), (A7) 



and the initial conditions A +1 v> = 1 and P+i,^ = 1- If we 
use Eq. (|A4[) the two additional coupled equations are 



— P-i M , = SNR V(r)A_ M , (A8) 

CtTf 

-^-A_i,^ = SNR iP(t)P-i,^, (A9) 

UTf 



with initial conditions A_i^ = —1 and P-i $ = 1. Note 
that since there is no relaxation, we do not need to 
have two equations for the ground state, we can just use 
d n P-i,^ = —SNR i/j(t)P-i^, but to keep the problem 
symmetrical we have decided to leave both equations in. 
This makes it easier to extend the theory to cases where 
upward jumps are possible. 



Thus to simulate Eqs. (|6.5p and (|6.6[) we simply solve 
the above two sets of two coupled differential equations 
and then combine them using Eqs. (|A1[) and (|A2[) . Note 
an equivalent derivation of these equations can be made 
by using the Kushner-Stratonovich equation [38| and 
then simply using Bayes theorem to invert these equa- 
tions for estimating unknown parameters [39l |. 
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